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S ■ Abstract 
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W). The appearance of high energy tau neutrinos due to — > i^r oscihations of 

^ ' extragalactic neutrinos can be observed by measuring the neutrino induced 

upward hadronic and electromagnetic showers and upward muons. We eval- 
uate quantitatively the tau neutrino regeneration in the Earth for a variety 
of extragalactic neutrino fluxes. Charged-current interactions of the upward 
tau neutrinos below and in the detector, and the subsequent tau decay create 
I muons or hadronic and electromagnetic showers. The background for these 

_ events are muon neutrino and electron neutrino charged-current and neutral- 

V/-^ ' current interactions, where in addition to extragalactic neutrinos, we consider 

■ atmospheric neutrinos. We find significant signal to background ratios for the 

. hadronic/electromagnetic showers with energies above 10 TeV to 100 TeV 

initiated by the extragalactic neutrinos. We show that the tau neutrinos from 
point sources also have the potential for discovery above a 1 TeV threshold. 
Q^' A kilometer-size neutrino telescope has a very good chance of detecting the 

appearance of tau neutrinos when both muon and hadronic/electromagnetic 
showers are detected. 
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H ■ I. INTRODUCTION 

A recent breakthrough in the study of neutrino oscillations came from the observation by 
the Super-Kamiokande experiment of a deficit of upward-going atmospheric muon neutrinos 
The observed electron neutrino flux was found to be consistent with the theoretical 
expectation from models of cosmic ray production of neutrinos. Furthermore, SuperK mea- 
surements are consistent with earlier experiments which detected anomalous ratios of 
the i>fj_ to Ue flux. The new high-statistics data disfavor scenarios in which i/^'s oscillate into 
sterile neutrinos (z/^) 0, and the data are consistent with z/^ to oscillation (99% CL) 
with a large mixing angle, sin^^ > 0.84 and a neutrino mass squared difference of 2 x 10~^ 
eV^ < Am^ < 6 X lO'^ eV^. 

Direct detection of Ur appearance is extremely difficult because at low energies, the 
charged-current cross section for producing a tau is small and the tau has a very short 
lifetime. Several long-baseline experiments with accelerator sources of z/^ f^HUl have been 
proposed with the goal of detecting tau neutrinos from oscillations, thus confirming the 
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SuperK results. The only convincing evidence of neutrino oscillations to date involves as- 
trophysical sources, neutrinos from the sun and atmospheric neutrinos. These observations 
involve indirect measurements, namely the disappearance of the expected neutrino fluxes. 

We have recently discussed the possibility of using a kilometer-size neutrino telescope 
to detect tau neutrinos from extragalactic sources of high-energy neutrinos such as Active 
Galactic Nuclei (AGN) and Gamma Ray Bursts (GRB), assuming z/^ <-> with the oscil- 
lation parameters of the SuperK experiment |T^. The probability for —>■ Vr is given by 

m 

In 2/l.27Am2(eV2)L(km)\ 
P{v, - L) = sm^2e sm^ ^(GeV) J ' 

Assuming two flavor oscillations, muon neutrinos produced in AGN or GRB would oscillate 
to tau neutrinos as they travel to the Earth. Over astronomical distances in the range of 
a megaparsec to thousands of megaparsecs, by measuring tau neutrino fluxes, one could, 
in principle, probe oscillations down to Am^ ~ 10"^'' eV^, nine orders of magnitude below 



current neutrino experiments ||l^Jl5|| . On the other hand, for the SuperK parameter range, 
with Am^ on the order of 10~^ eV^ and sin^ 26' ~ 1, the oscillation probability is 0.5. It is 
this latter possibility that we explore in this paper. 

We use the simplest assumption for the flavor content of extragalactic sources of neutri- 
nos, in the absence of oscillations, for the ratio z/g : z/^ : z/,- to be 1 : 2 : 0. This is based 
on a counting argument applied to vr — /iz/^ and fi — > t'^ez/g processes. With the two-flavor 
oscillations suggested by the SuperK experiment, the flavor ratio becomes 1:1:1 after the 
neutrinos have traveled over astronomical distances. Even in the three-flavor oscillation sce- 
nario, the ratio is still 1:1:1, because the path length of high energy extragalactic neutrinos 
is much larger than any neutrino oscillation length supported by the solar, atmospheric or 



accelerator data [p!6|. 



The ratio for z/g : z/^ might get modified at high energies due to muon cooling [0 . In ad- 
dition, z/g from neutron decay might give significant contribution, resulting in neutrino fluxes 
dominated by electron neutrinos as in the case of diffuse neutrino fluxes from propagating 



cosmic rays |18[. We comment qualitatively in the discussion section on how our results 
are altered with more realistic, flavor- dependent neutrino energy cutoffs. Regardless of the 
flavor content of the source, the maximal mixing suggested by the SuperK experimental 
results mean that there will be an appreciable tau neutrino component at the Earth, so one 
is interested in tau neutrino detection in high energy neutrino telescopes such as ANTARES 



TP| , NESTOR [^, AMANDA |^| and the next generation of large underground detectors 
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Tau neutrino detection requires an understanding of the effect of propagating through the 
Earth on the tau neutrino flux. The propagation of ultra-high energy tau neutrinos through 
the Earth is quite different from muon and electron neutrinos. The Earth never becomes 
opaque to tau neutrinos, while muon and electron neutrinos are absorbed via charged- 



current interactions before reaching the opposite surface [|14[. Ultrahigh-energy tau neutrinos 
interact in the Earth producing taus which, due to the short lifetime, decay back into 
tau neutrinos with lower energy. This cascade continues until the tau neutrinos reach the 
detector on the opposite side of the Earth or until the energy of the neutrinos is small enough 
that the interaction length of the neutrino is longer than the path length through the Earth. 
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The energy and nadir angle dependence of the extragalactic tau neutrinos fluxes have been 
examined quantitatively in Refs. |12,^|. For certain fluxes, those that do not decrease too 



steeply with energy, there are significant enhancements of the tau neutrino flux relative to 
the muon neutrino flux at energies below ~ 10^ GeV. 

An enhancement of the tau neutrino flux does not necessarily translate to dramatic 
modifications of the standard model (no-oscillation) rates for upward-going muons, espe- 
cially in view of uncertainties in the normalization of the extragalactic fluxes. However, by 
comparing rates for upward-going muons with rates for upward hadronic/electromagnetic 
(EM) showers, the signature of tau neutrino interactions is unambiguous for a large range 
of neutrino fluxes. In the next section, we briefly introduce the extragalactic and generic 
Vt- fluxes ~ i?"" for n = 1, 2 that are used here. After reviewing neutrino propagation 
through the Earth, we describe Vr signatures. The fluxes considered here have a range of 
energy behaviors. Even if the normalizations of the neutrino fluxes are uncertain, and in 
some cases optimistic, it is useful to make quantitative comparisons of the event rates for 
upward muons and upward hadronic/EM showers, with and without neutrino oscillations, 
which we do in Section IV. The quantitative results for specific models lead to model in- 
dependent conclusions, which we summarize graphically. Tau neutrino appearance would 
provide an independent confirmation of the SuperK results and would point towards the 
better understanding of physics beyond the Standard Model. 



II. HIGH ENERGY NEUTRINOS SOURCES 

Active Galactic Nuclei are the most luminous objects in the Universe. Most of this 
radiation comes from their central region, indicating that the energy radiation most likely 
comes from accretion of matter into a superheavy black hole. Protons within the AGN 
may get accelerated via first order Fermi acceleration to very high energies. They interact 
with protons and photons in the infalling gas, or they may exist in the jets along the 
rotation axis and interact with photons there. Photon-proton and proton-proton interactions 
produce pions, which decay into charged leptons, neutrinos and photons. Energetic photons 
{E^ ~ 100 MeV) from about 40 AGN observed by the EGRET collaboration 0, and TeV 
photons have been detected from Mkn 421, Mkn 501 fgD and lES2344-h514 pg. Although 



these photons are conventionally explained by inverse Compton scattering from energetic 
electrons, this explanation is not without problems, and a hadronic origin of gamma-ray 
photons from AGN is a viable alternative [27]. If a large fraction of the observed energy 
in high energy photons from AGN is produced in hadronic interactions, then AGN are also 
powerful sources of ultrahigh-energy (UHE) neutrinos [P^-pi . 



In Fig. 1 we show neutrino fluxes predicted in the AGN models of Stecker and Salamon 
28| and Mannheim Model A [29|. Both of these models predict neutrinos fluxes that rep- 



resent the upper bounds for their class of the models. In particular, the Stecker-Salamon 
flux is an upper bound for AGN core emission, while Mannheim Model A is an upper bound 
for AGN jet emission models. Stecker-Salamon flux is bound by the the diffuse X-ray back- 
ground, while Mannheim flux is bound by the extragalactic gamma ray background. The 
steep, low energy neutrino flux in Mannheim's model is the emission from the host galaxy 
via pp interactions of the AGN protons in the galactic gas disk. Since this part of the flux 
is derived with the assumption that all protons end up in the disk, it should be regarded 
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as an upper bound. Stecker-Salamon flux at energies above 1 PeV may get reduced due to 
the cooling of pions and muons in the strong magnetic flelds of AGN cores [|l7l. The fluxes 
plotted in Fig. 1 are for the sum of muon neutrino plus antineutrino, at the source, namely, 
without accounting for oscillation over astronomical distances. We label the fluxes in the 
absence of oscillations by -F,^+p. 

Another extragalactic source with powerful radiation and possibly associated high energy 
neutrino flux are the gamma ray bursts (GRB). Several models have been proposed in order 
to explain the origin of GRB's ||3^j33[| . In the flreball model |^3|, the gamma ray bursts are 



produced by the dissipation of the kinetic energy of the relativistic expanding flreball with a 
large fraction, > 10%, of the flreball energy being converted by photopion production to high 



energy neutrinos ||3^. Photomeson production takes place when extremely energetic protons 
accelerated at high energies in the ultra-relativistic shocks interact with synchrotron photons 
inside the flreball. The decay of these charged pions and subsequently produced muons then 
produce electron and muon neutrinos. Contributions from proton-proton collision can be 
neglected in this model. In Fig. 1 we show the neutrino fluxes for gamma ray burst model 



of Waxman and Bahcall (GRB_WB) I^Sj, in which they parameterize the flux by 

= 4.0 X lO""^-'^, 

where a = 13 and = 1 for £^ < 10^ GeV and a = 8, = 2 for 10^ < £" < 10''' GeV and 
a = 1, n = 3 for E > 10^ GeV. 

Theoretical work has been done to set upper bounds on high energy neutrino fluxes from 
AGN jets and GRB |^. The bounds are based on the theoretical correlations between 
the cosmic ray flux and/or the extragalactic gamma ray flux and the neutrino flux. These 
bounds have some model dependence, and they tend to be weaker in the range of energies 
considered here than at higher energies {E ~ 10^ — 10^ GeV). The AGN and GRB neutrino 
fluxes used here satisfy these bounds. 

Cosmic topological defects (TD) such as magnetic monopoles, cosmic strings and domain 
walls are predicted to be formed in the Early Universe as a result of symmetry breaking 
and phase transition in Grand Unifled Theories (GUTs) of particle interactions. In the 
TD models, 7-rays, electrons (positrons), and neutrinos are produced directly at ultra- 
high energies via cascades initiated by the decay of a supermassive elementary "X" particle 
associated with some Grand Unifled Theory, rather than being produced in high energy 
hadronic interactions. The X particle is usually thought to be released from topological 
monopoles left over from GUT phase transition. It decays into quarks, gluons, leptons. In 
this paper, we consider neutrino fluxes from topological defects models of Sigl-Lee-Schramm- 
Coppi (TD_SLSC) |6[ and the model of Wichoski-MacGibbon-Brandenberger (TD_WMB) 
P7| . The main difference between these two models is the main channel for energy loss of 
the string network, in the former it is the gravitational radiation, while in the later it is the 
particle production. Both of these fluxes should be regarded as upper limits for TD models, 
because they have been constructed in such a way to satisfy the bound imposed by the 



measured cosmic ray and gamma ray fluxes |38|. These fluxes are shown in Fig. 1, where we 



take representative flux of WMB model with the string mass parameter giving the largest 



neutrino flux that is consistent with cosmic ray data. This flux is also below the Frejus ||39 
and Fly's Eye experimental limits on the neutrino flux. 



We also consider two generic fluxes that have a power law behavior. The flux 
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(2) 



gives numerically stable results, however, our calculations with a flux with 
unstable at very high energies. Consequently, we use 



F:^,{E) = 10-"(i?/GeV)-i 



1 



(1 + E/108GeV)2 



fcm-^s-^sr-^GeV-^l 



E'^ is 



(3) 



as a way to cutoff the high energy behavior. We show results for attenuated fluxes for 
neutrino energies up to 10^ GeV. We have chosen the multiplicative factors in E~^ and E'"^ 
fluxes in such a way that they exceed atmospheric flux at neutrino energies between 10 TeV 
and 100 TeV. The upper bound for strong source evolution discussed recently by Waxman 
and Bahcall [^] would correspond to a limit of 2 x lO^'^-E'"^ (in the same units as Eq. (2)), 
a factor of 5 smaller than the choice of normalization we take in this paper. 

We also show in Fig. 1 the atmospheric neutrino flux at zenith angle of 0° and the 
horizontal flux. In our evaluation of the atmospheric backgrounds, we use the atmospheric 



muon and electron neutrino fluxes as a function of zenith angle [41 



We do not consider the neutrino flux from cosmic ray interactions with the microwave 
background. This diffuse neutrino flux is typically present at energies higher than we consider 
here and it gives low event rates [^. Furthermore, the cosmic ray interactions with 
the solar atmosphere are another source of neutrinos, however, for energies above a TeV, 
the flux scales as E~^ or steeper As we see below, tau neutrino regeneration will not 



be a very important feature in fluxes with such large spectral indices. 



III. TAU NEUTRINO PROPAGATION THROUGH THE EARTH 

For neutrino energies above 1 TeV, the oscillation probability for i/^ — > Ur in the Earth 
is less than a percent for the parameters constrained by the SuperK experiment. As a con- 
sequence, we can neglect neutrino oscillation in our evaluation of tau neutrino propagation 
accounting for interactions in the Earth. 

The coupled transport equations for the fluxes of the tau neutrino and its charged partner 
are given by 



dF,^{E,X) F,^{E,X) 



dX K^{E) 

roo 

+ dEyGr^Ur {Ey, E, X) 

and for taus as, 

dF,{E,X) F,{E,X) F,{E,X) 



/•oo 

/ dEyG,^^,^{Ey,E,X) (4) 

J E 



dX K{E) \f^{E,X,e) 

roo roo 

+ / dEyGr-^r{Ey, E, X) + / dEyGy^^r{Ey,E^X). 
J E J E 



(5) 



Here Fy^{E,X) and Fr{E,X) are differential energy spectrum of tau neutrinos and tau 
respectively, for lepton energy E, at a column depth X in the medium deflned by 
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X= p{L')dL'. (6) 
Jo 

The density of the medium a distance L from the Earth-atmosphere boundary, measured 
along the neutrino beam path, is p{L). The lepton interaction length (in g/cm^) is X{E) 
and Xf'^{E, X, 9) is the decay length of the tau. 

The functions Gi^j schematically represent interaction or decay contributions to lepton 
j from lepton i. We limit our evaluations of the tau neutrino flux to E^^ < 10^ GeV. 
Consequently, we can ignore several terms in the coupled differential equations: the term 
with Gr^T and the term —Et/Xt, both in Eq. (5). Only the decay contribution to the last 
term in Eq. (2) (Gr^^^) is included in our evaluation. This is justified by the fact that the 
tau decay is significantly more important than interactions for the energy range of interest, 
namely E < 10^ GeV. The neglected terms start contributing for lepton energies on the 



order of 10® GeV. Detailed formulae for Gi^j appear in Ref. |T2 



In our previous work, we have described an analytic method for solving these transport 
equations ll2|, based on the method of Naumov and Perrone |^4|. We have evaluated the 



upward Ur flux for a selection of initial fluxes |jT2|. We have shown that for "fiat" initial 
neutrino fluxes {E ~ E~^), a significant number of high energy i/^'s cascade down in energy, 
resulting in enhanced low energy flux relative to the attenuated z/^ flux. Here, we evaluate 
the tau neutrino flux for a more comprehensive selection of incident fluxes, including both 
neutrino and antineutrino attenuation. In all of our results below, we evaluate the sum of 
neutrino plus antineutrino fluxes or rates. 

In Fig. 2 and 3, we show the attenuated tau neutrino plus antineutrino flux (blue line) 
and attenuated muon neutrino plus antineutrino flux (Red line), scaled by a factor of the 
neutrino energy E, assuming the equal fluxes of tau neutrinos and muon neutrinos incident 
on the surface of the Earth at a nadir angle of 0° for E^^^{Ey) = 0.5 x 10~^^E~^, F°^_p(-Ey) = 
0.5 X 10"^^-^ the Stecker-Salamon AGN model, the Mannheim AGN (Model A), the two 
topological defects models, the Waxman-Bahcall GRB model and the atmospheric flux. We 
note that the enhancement of the tau neutrino flux relative to the initial flux and also 
to the muon neutrino flux, is prominent for the flat fluxes, such as E°{E^) ~ E~'^, the 
Stecker-Salamon AGN model and the topological defects model of Sigl et al. In case of the 
atmospheric flux, which represents the background, the enhancement is very small due to 
the steepness of the initial neutrino flux. 

The angular dependence of the upward flux is also distinct. As an example, in Fig. 
4, for the AGN model of Mannheim (Model A) |2^ , we show the ratio of the neutrino flux 
scaled by the flux at X = for two nadir angles, ^ = 0° and 9 = 30°, as a function of 
neutrino energy. Because of the shape of the initial flux, steep for energies below 10^ GeV, 
and flat for higher energies, the enhancement of the tau neutrino flux becomes significant 
only for energies above 10® GeV. At fixed energies of 10^, 10^, and 10® GeV, we show the 
same flux ratios function of nadir angle. 



IV. DETECTION OF i/^ APPEARANCE 

Detection of muon neutrinos, in general, is via their charged-current interactions. Pro- 
duced muons have very large average range making the effective volume of an underground 
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detector significantly larger than the instrumented volume. On the other hand, tau neu- 
trino charged- current interactions produce tau, which has a very short lifetime, making its 
detection extremelly difficult. Only at very high energies, > PeV, the production and 
decay vertices are separated by a measurable distance providing a distinctive signature of 
tau neutrinos ("double-bang" events) However, the predicted neutrino fluxes are low 
at these energies. For z/t-'s in the energy range of 10^ — 10^ GeV considered here, the pro- 
duced tau decays after a very short pathlength back to Ur plus leptons or hadrons. Tau 
neutrinos will interact via neutral currents, producing a hadronic signal as well. Therefore, 
the signals of tau neutrino interactions below the double-bang threshold are muons from tau 
decay, or hadronic/EM showers from the tau production and/or decay. In the first case the 
background to tau production of high energy muons is z/^ charged-current interactions. In 
the latter case, the backgrounds are z/^ neutral current and Ue charged-current and neutral 
current interactions. The background rates shown below are with the assumption that the 
electromagnetic shower from ^ e charged current interactions cannot be distinguished 
from hadronic showers. As a consequence, we evaluate the hadronic/electromagnetic (EM) 
shower rates. 

We assume in the analysis presented below that the z/^ charged-current events and I'r — ^ 
T fi events can be rejected from the contained hadronic/electromagnetic shower signal. 
In both cases there is a hadronic shower which includes muons, however, the muons in the 
hadronic showers from pion and kaon decays are significantly less energetic than the muons 
from z/^ ^ /i and Ur t ^ fi. In the latter case, the energy of the shower is ~ 1/2 
the incident neutrino energy E^, and the energy of the muon is ~ 1/6 — 1/2 E^. On the 
other hand, muons coming from particle decays in the hadronic shower are considerably less 
energetic because of large particle multiplicities. The average charged particle multiplicities 
for hadronic interactions at -^i > 40 GeV are larger than ~ 10 particles so individual 
muon energies from charged pion and kaon decays are less than ~ 5% of the incident neutrino 
energy. The hadronic shower and very energetic muon of the "muon signal" should stand out 
in comparison to the hadronic/EM shower signal in a detector with good energy resolution 
like the proposed kilometer-cubed detector IceCube 

We describe the evaluation of the muon and hadronic/EM shower event rates. The event 
rates for z/^ r ^ and z/^ — are evaluated and compared with the no oscillation rates. 
We evaluate the hadronic/EM shower rates for signal and background, then compare with 
the hadronic/EM shower rates assuming no oscillations of z/^. The relative rates of muons 
and hadronic/EM showers prove to be the most effective diagnostic to neutrino oscillations 
with the SuperK parameters. 



A. Muon Event Rates 

The standard evaluation of the muon event rate per solid angle for neutrino interactions 
with isoscalar nucleons N {v^N — > ^X) follows from the formula |^ 



Rate= ANa H dE, f dy{R,{E,{l - y), i^;--)^^^^-^ (7) 
j£;min J ^ ay 



xF,{E,,X)Q{E,{l-y)-E^ 
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where y is the neutrino energy loss, y = {Ey — E^)/Ey, and dacc{Ey,y) / dy is the charged 
current differential cross section. Fy{Ei,,X) is the upward neutrino or antineutrino ffux 
which depends on angle implicitly through the pathlength X. We assume that the initial 
ffuxes of muon neutrinos and antineutrinos that reach the Earth are equal, their sum in the 
oscillation scenario being half of the muon neutrino plus antineutrino flux produced at the 
source. The fluxes of neutrinos and antineutrinos at the detector are different because of the 
difference in charged and neutral current cross sections below energies of 10^ GeV PD|,ET 



however, for these energies and fluxes, the antineutrino event rates differ for the neutrino 
event rates by at most ±20%. The average range of a muon, is the range 

of a muon produced in a charged- current interaction with energy E^ which, as it passes 
through the medium, looses its energy via bremsstrahlung, ionization, pair production and 
photonuclear interaction and arrives in a detector with an energy above -E"™™. Avogadro's 
number is and A is the effective area of the detector. All of the event rates calculated 
are for the sum of neutrino plus antineutrino contributions to /i^ + /i~ production. 

The rate for muons produced by the tau neutrino charged current interactions followed 
by the tau leptonic decays is given by a modified equation, taking into account the branching 
fraction for r Vt^^^^ and the decay distribution of the muon via dn{Er)/dz, where z = 
E^/Et-. The decay formulae used here are listed in the Appendix A. The differential event 



rate is 



Rate= ANa j^,^ dE^ j dy j dz{R^{Ey{l - y)z, j^min^^ ^^(-^^0-^ v)^) (g) 



dcTccjEu, y) 
dy 



X ^^^J^Fy{Ey, X)e{Ey{l - V) Z > E^^) 



In Fig. 5 we show the neutrino processes that contribute to the muon production. In our 
evaluation of the muon event rates we use the Earth densities of the Preliminary Earth Model 
(PREM) described in Ref. [|^. We have used the PREM to determine an average density 
for a given nadir angle, then used that average density to evaluate the attenuated fluxes. 
We use the muon range evaluated by Lipari and Stanev |^9|. The neutrino and antineutrino 



cross sections have been evaluated using the CTEQ5 parton distribution functions [Q. The 
effective area A is taken to be 1 km^. 

In Figs. 6-9, we show muon event rates for ~ E^^, F°^_p ~ E"^, AGN_SS, 

AGN_M95, TD_WMB, TD_SLSC and GRB_WB for Ef"" = 1, 10, 100 TeV. Blue lines 
correspond to the upward + fi~ events from z/,- + z/,- + z/^ + charged-current interactions 
(including t ^ fi decay), while the red lines are the background contribution from z/^ + z/^ 
charged-current interaction only. We note that the muon enhancement due to the tau neu- 
trino contribution for E~^ flux is almost factor of 2 for small angles and 25% for large angles 
with E"™™ = 1 TeV. The enhancement is less pronounced at small nadir angles for increasing 
threshold energies, for example, the blue line is about 60% enhanced relative to the red line 
for -E™" = 10 TeV for the E~^ flux in Fig. 6. A similar enhancement is present for the 
TD_SLSC. For steeper fluxes, such as AGN_SS, AGNJVI95 and the enhancement due 
to tau neutrino contribution is much smaller, of the order of 20-25%. 

The muon event rates from the atmospheric neutrino background are shown in Fig. 



10a). The input flux is the angle dependent muon neutrino flux of Agrawal et al. 



The atmospheric tau neutrino flux is very low, as the tau neutrinos are produced in the 
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atmosphere by comic ray interactions with nuclei in the atmosphere, which produce Ds 



whose leptonic decay, Ds tv^, gives v-r I^T]. The rates for the atmospheric tau neutrinos 



are shown in Fig. 10b). In the evaluation of the event rates, we neglect oscillations of 
atmospheric neutrinos as they travel to the Earth and the oscillations through the Earth 
since the oscillation probabilities are small above our minimum energy of 1 TeV. 

The atmospheric neutrino events represent a background for detection of extragalactic 
neutrinos. For a muon energy threshold of 1 TeV, the background is large, 400 — 2000 events 
per year per steradian for 1 km^ effective area detector. For a muon threshold of -E™™ = 10 
TeV, the event rates range between 6 — 80 events per year per steradian. A comparison of 
the event rates from Figs. 6-9 with the atmospheric muon neutrino background indicates 
that detection of neutrinos from AGN might be possible with E^^'^ = 10 TeV or 100 TeV. 

The rates for muon events shown in Figs. 6-9 come from assuming that the tau neutrino 
and muon neutrino fluxes are equal and are half the flux of muon neutrinos produced at 
the source. Testing the oscillation hypothesis with muon neutrinos alone will be difficult. 
We see that with the exception of the and TD_SLSC fluxes, the observed muon rate is 
about half of what one would expect in the absence of oscillations. Given the uncertainties 
in the normalizations of the predicted fluxes, this factor would not unambiguously signal the 
presence of tau neutrinos from oscillations. The situation with the and TD_SLSC fluxes 
is only slightly better. There, in the oscillation scenario, the measured muon event rate is 
about 80% of the no oscillation prediction at = 0°, but less than 70% of the prediction for 
horizontal events. Testing the oscillation hypothesis by measuring upward muons only will 
be very difficult. 

The relatively small contribution to the muon rate from z/^-'s, despite the fact that the 
attenuated flux of tau neutrinos is larger than that of the muon neutrinos, is due to the fact 
that the muon carries a small fraction of the initial tau neutrino energy. Consequently, for 
a muon of a given energy, if it comes from a tau neutrino (which interacted producing a tau 
that subsequently decayed to a muon), the initial tau neutrino has a much higher energy 
than a muon neutrino which produces a muon directly via the charged current u^^N —>■ fiX 
process. All predicted neutrino fluxes decrease with energy. Even with some "pileup," the 
tau neutrino fluxes are decreasing fast enough that the muon energy fraction results in 
sampling a much smaller tau neutrino flux than the corresponding muon neutrino flux. It is 
this observation that leads one to consider signals that carry a much larger fraction of the 
incident tau neutrino energy. 



B. Upward Hadronic/Electromagnetic Showers and Their Detection 

The hadronic/EM shower signal of Ur interactions is a much more promising flnal state 
from the theoretical point of view than the muon signal. The beneflt is that the hadronic 
showers include both production hadrons and tau decay hadrons, so there is a much higher 
fraction of the incident tau neutrino energy visible in the detector |^ . The next generation 



of neutrino telescopes may not be able to distinguish between hadronic and electromagnetic 
showers, so we include in the signal and in the background, processes that include hadrons 
and electron. As mentioned above, we assume that the high energy muon associated with 
the target jet in charged current interactions will be used to veto the process u^^N — > fiX. 
Distinguishing electromagnetic from hadronic showers might be possible by looking at the 
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difference between the front to back ratio of tlie cascade Cherenkov light, and perhaps by 
the number of residual vr ^ /i ^ e decay, although this is considered to be experimentally 



difficult 153 



The processes that go into our evaluation of — >• hadrons are 

Vt-N ^ T + hadrons, t ^ + hadrons , 
PrN r + hadrons, r ^ z/,- + e + i^e , 
v^N Vt- + hadrons . 

For the charged-current interactions, the hadronic/electromagnetic energy is the sum of the 
energy carried by the hadrons in tau production, as well as the tau decay hadronic energy 
or tau decay electron energy. 

The background for the hadronic/electromagnetic showers is due to the z/^ and z/g neutral 
current interactions, and Ve charged-current interactions are 

^'M.e + N ^ u^^e + hadrons , 
i>e + N ^ e + hadrons . 

For the Uf, flux, we assume it is equal to the flux in the SuperK oscillation scenario. All 
of the processes that contribute to the hadronic/EM showers are shown in Fig. 11. 

The tau neutrino shower event rate per unit solid angle from charged-current interactions 
followed by the tau hadronic decay is given by 

Rate= VN,£jKldyldz'^'^^^^^^f^K,(E.,,X) (9) 



xe(E,Av + {1 ~ v){i - z)) - EZ") 



The hadronic energy from the broken nucleon i^g^* = E,yy and the hadronic energy from 
the decay E^^^^^ = Ei^{l — y){l — z) are added to get the total shower energy. Again, 
y = {El, — Et)/ Ey for incident neutrino energy Ey, while z = E^/Er, where E'^ is the energy 
of the neutrino from the tau decay. The differential distributions for the hadronic decay 
modes are shown in the Appendix A. For the electronic decay of the tau, the differential 
distribution dn/dz is replaced by the purely leptonic distribution in terms of z' = E^/Et-. 
The theta function is replaced by 

e{Ey^{y + (1 - y){l - z)) - Etn - ^iKAv + (1 - v)^') " ^sT) • (10) 
The neutral current background event rate is given by 

Rate = VNa H dE, f X)e{Eyy - E^^ , (H) 

J dy 



shr 



while the electron neutrino charged current background rate is given by 

Rate = VNa H dE, [ dy^^^^^^^Fy{E,, X)Q{E, - E^"). (12) 
JET" J dy 



dy 

In Figs. 12-15 we show the upward hadronic/EM shower event rates as a function of 
the nadir angle for Es^r > E^^ where E^^ = 1 TeV, 10 TeV and 100 TeV for input fluxes: 
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FO+p ~ E-\ F^+p ~ AGN_SS, AGN_M95, TD_WMB, TD_SLSC and GRB_WB, all 
assuming that V = 1 km^. The blue lines correspond to the event rates from z/^ + z/^ + z/g + z/g 
charged- current interactions (and r — >• i't-+ hadrons decay) and from i/T- + z/T- + i/^-|-i/^-|-z/e + i^e 
neutral current interactions. The red lines are the contributions from i/^ + 1/^ neutral current 
interaction and Ue + Pe charged and neutral current interactions. We do not include + 
charged-current interactions in our calculation because these events can be vetoed by the high 
energy muons produced in the interactions. All of the rates shown in these figures assume 
equal neutrino and antineutrino fluxes. They are performed in the oscillation scenario where 
the ratios of the fluxes '■ '■ are 1:1:1. 

Prom Fig. 12a) we note that in the case of the E'^ flux, the contributions from tau 
neutrinos are large, a factor of 4 times larger than the muon neutrino plus electron neutrino 
contribution at zero nadir angle. For horizontal showers, the enhancement factor is smaller, 
about 2 for all the energy thresholds that we consider. Similarly, for flux, the tau 
neutrino contribution is a factor of 1.7 times larger than the muon neutrino plus electron 
neutrino contributions for upward showers. 

Similar conclusions can be drawn from the plots of the other fluxes. The shower event 
rates including z/^ + z/^ + z/^ + z7^ + z/e + z/e are significantly enhanced relative to the rates 
from z/^ + z/^ + z/e + z/g in the oscillation scenario. The AGN_SS rates at zero nadir angle 
are comprised of 60% tau neutrino induced, decreasing to about 40% tau neutrino induced 
for horizontal showers, as shown in Fig. 13a). AGN_SS flux gives 25-80 shower events for 
= 10 TeV and 6-45 events for E^^ = 100 TcV with negligible atmospheric background. 
In Fig. 13b) we show event rates for AGN_M95 model. We find 3-6 shower events per year 
per steradian for = 10 TeV, with atmospheric background of 2-16 events. Detection 

of events with higher energy threshold would require looking at almost horizontal showers, 
where the background is small. 

The TD_WMB model in Fig. 14a) shows an enhancement of between 2.1-2.3 for zero 
nadir angle, and a factor of 1.7 for almost horizontal showers. Fig. 14b) shows the more 
striking enhancement in the TD_SLSC model, where the enhancement is a factor of between 
3.7 to 6.2 at zero nadir angle, to a factor of 2 for large nadir angles. However, due to the 
particularly low normalization of the TD_SLSC flux, the kilometer-size detector would not 
be sufficient for its detection. Fig. 15 shows the GRB_WB model in which the enhancement 
factor is between 1.5 to 2, depending on energy threshold and angle. The event rates for 
showers with energies above 10 TeV are comparable with the background, but higher energy 
threshold of 100 TeV would still give a few events per year for large nadir angle with negligible 
background. 

Fig. 16a) and 16b) show the shower event rates for the atmospheric z/^ + z/^ + i^e + Pe 
and Vr + fluxes, respectively. For showers with energies above 10 TeV, the event rates 
are twice as large as for the E^^ flux at small nadir angle. For E'^^^ = 10 TeV, we flnd the 
event rates for the showers to be about 8-18 per km^ per year per steradian for the E^^ flux, 
compared with the atmospheric background of 2-16. 

The AGN rates will stand out above the atmospheric background for -E™" ~ 10 TeV. The 
GRB_WB rates are more than half of the atmospheric neutrino rates at the 10 TeV shower 
threshold at small nadir angles. The TD rates are all quite low overall, and in comparison 
to the atmospheric background rates. 

Since one does not measure separately the tau neutrino induced shower rates and the 



11 



nmon and electron induced shower rates, given a particular model, one can compare the 
rates with the oscillation hypothesis to the predicted rates without oscillations. We discuss 
the rates for specific models here, then discuss a more model independent analysis in the 
next section. To illustrate the effect of oscillations we plot in Figs. 17-20 the ratio of the 
shower event rates from i/,-, and Ve in the oscillation scenario to the shower rates from 
and z/g in the standard model, with no oscillations. 

From Fig. 17a) we note that for the flux, the shower event rates in the oscillation 
scenario are a factor of 3.3-3.7 larger than in the no oscillation case for i^'shr" = 1 — 100 TeV 
for 6 — 0°. They are a factor of 1.6 enhanced for the horizontal shower rate. For the E~'^ flux, 
the enhancement is a factor of 1.4-1.6 relative to the no oscillation case for — 1 — 100 
TeV, shown in Fig. 17b). In the case of AGN models, if one assumes oscillations, the shower 
event rates are factor of 1.8-2.1 larger for AGN_SS at zero nadir angle, decreasing to 1.5 for 
nearly horizontal showers, as shown in Fig. 18a). Fig. 18b) shows a ratio ranging between 
1.4-1.9 for AGN_M95 for small nadir angles. 

From Fig. 19a), we note that the shower event rates for TD_WMB are factor of 1.8-2.1 
enhanced for energy thresholds of 1-100 TeV for the upward neutrinos, while the enhance- 
ment is a factor of 1.5 for almost horizontal showers. In the TD_SLSC model. Fig. 19b), 
the shower event rate is a factor of 3-3.6 enhanced at small nadir angles, and a factor of 
1.6 enhanced for horizontal showers. In the case of GRB_WB model. Fig. 20 shows an 
enhancement of 1.4-1.7 for ET^ = 1-100 TeV. 



C. Relative Rates 



We have shown that comparison of the muon and shower rates serves as a diagnostic for 
Ut oscillations over astronomical distances. For example, for the E'"^ flux. 



Ratio 



{vr + i'n + iye^ shower)osc/ (i^ix + i^e^ shower)r 



~ 1.5 



while 



Ratio 



~ 0.5 . 



(13) 



(14) 



The ratios include contributions to showers and muons from antineutrinos. This feature of 
a deflcit of muon rates and an excess of shower rates in the oscillation scenario compared 
to the no-oscillation scenario is a generic feature of all the neutrino spectra in Fig. 1. To 
demonstrate this point quantitatively, we deflne a ratio of ratios. 



R = 



Ratio 


{ur + i^ij. + i^e^ shower)osc/ (i^ij, + i^e^ shower)no-osc 


Ratio 


(^T ~l~ ^fj, 1^'^ osc / iy fj, A'')no— osc 





(15) 



As Figs. 6-9 and 17-20 illustrate for individual fluxes, R depends on energy threshold 
and angle. We show in Fig. 21 the band of R spanned by the representative models of Fig. 
1 for three thresholds in muon or shower energy: a) 1 TeV, b) 10 TeV and c) 100 TeV. We 
note that R depends on nadir angle and threshold energy, however, it! ^ 2.4 independent of 
the initial flux. For a given model, measured rates will be very distinct from predicted rates 
if the SuperK results for oscillation parameters are correct. 
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Determining 



Ratio 


(shower ratc)i„easured/(/^ rate)^^^^^^^j 




Ratio 


{l^H + l^e^ shower)„o_osc/(l^M A*)no-osc 



(16) 



nevertheless rehes on theoretical input for the "no-oscillation" flux. Different energy behav- 
iors of incident fluxes will have implications for the angular and energy dependence of the 
event rates of upward muons and upward hadronic/EM showers, allowing for an indirect 
characterization of the energy dependence of the source. 

A more model independent test of the oscillation scenario would be to compare the 
measured ratio of showers to muons with the no-oscillation predictions, on an absolute 
scale. This requires a crude separation of the different energy behaviors of the fluxes of Fig. 
1. By only looking at the ratio of shower to muon events, one could confuse the AGN_M95 
no-oscillation ratio with the similar GRB_WB oscillation ratio. However, experimentally, 
the energy and angular dependence of the muon event rates for the two fluxes are quite 
different, and GRB neutrinos would reveal themselves by time correlations to observed 
ORB events. One category of fluxes, with not too steep energy behavior {E~^, AGN_SS, 
TD_WMB, TD_SLSC and GRB_WB), have a reduction in the event rates from a muon 
threshold of 1 TeV to 10 TeV at nadir angle of 0° by a factor of less than 4, whereas for the 
other two steeper fluxes {E~^ and AGN_M95), the reduction is by more than a factor of 6. 
These reduction factors are independent of whether or not oscillations occur. The steeper 
fluxes are also distinguished by muon event rates with a less marked dependence on nadir 
angle. If one separates the steep from the less steep examples used here, then at all three 
threshold energies, the band of oscillation shower to muon ratios does not overlap with the 
band of no-oscillation ratios. This is shown graphically in Figs. 22 (a-f) for 1 TeV, 10 TeV 
and 100 TeV thresholds, respectively. In fact, for the 100 TeV threshold, one does not need 
any information about the energy dependence of the initial flux, but in this case, the event 
rates are expected to be low. 



V. DISCUSSION 

We have studied signals for — > i/^ oscillations with extragalactic high energy muon 
neutrinos. Assuming SuperK oscillation parameters, muon neutrinos convert into tau neu- 
trinos as they travel megaparsec distances, with both fluxes being equal at the surface of the 
Earth. High energy muon neutrinos get absorbed as they pass through the Earth, while tau 
neutrinos cascade down to lower energies. We find this enhancement of the fiux in the 
low energy region to be prominent for flat initial spectrum, such as E~^, the AGN model 
of Stecker and Salamon, and the topological model of Sigl et al. For steeper spectra, the 
enhancement is small because the number of higher energy neutrinos that contributes to 
the lower energy flux via tau decay is relatively small compared to the low energy flux of 
neutrinos. 

Upward tau neutrinos, once they reach the detector, interact producing tau Icptons 
which decay with very short lifetimes. We have considered muons from tau decay as well 
as its hadronic decay mode. Since the planned detectors are unable to distinguish between 
hadronic and electromagnetic showers, we have included all the processes that give both 
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hadronic and electromagnetic showers. We find that upward muons alone would not be 
sufficient to separate the tau neutrinos contribution, due to the large background from z/^ 
charged- current interactions, the small branching fraction for r — >■ /i decay mode and the 
model uncertainty for the incident neutrino flux. 

In the case of upward hadronic/EM showers, we find that tau neutrinos give significant 
contributions, signaling the Vr appearance. Given the uncertainties in the normalizations 
of the extragalactic neutrino fluxes, combining muon rates and hadronic/EM rates offer the 
best chance to test the — >■ Vj. oscillation hypothesis. 

As concluded in earlier work [^,TT|, in general, an energy threshold of between 10 TeV 
and 100 TeV for upward muons and showers is needed in order to reduce the background 
from atmospheric neutrinos. We find that diffuse AGN neutrino fluxes, as described by the 
Stecker-Salamon and Mannheim models, as well as neutrinos from GRBs can be used to 
detect tau appearance. By measuring upward showers with energy threshold of 10 TeV, and 
upward muons, the event rates exceed the atmospheric background and are about a factor 
of 1.5-2 larger than in the no-oscillation scenario. 

Here we also comment on the effect of muon and pion cooling to the flavor ratio. Athar et 
al. in Ref. [0 have shown that with a negligible electron neutrino content at the source, the 
electron neutrino content at the Earth (in the three-flavor model) is reduced if not negligible 
compared to the nearly equal muon and tau neutrino fluxes. Keeping the energy spectrum 
unchanged, this means that the hadronic/electromagnetic shower background, which has 
significant contributions from v^N eX with Ue > would be reduced. Electron (anti- 
) neutrinos from processes in the propagation of cosmic rays may dominate at some energies 
||r8| . We have not considered that possibility here because of the low rates below 1 PeV. 

Steepening of the energy spectra displayed in Fig. 1 due to a neutrino energy cutoff from 
pion and muon cooling will have implications for the tau neutrino 'pileup', especially for the 
flatter spectra where the pileup is more pronounced. As an estimate of the lower bound on 
the relative enhancement of the hadronic/EM signal compared to the muon signal, one can 
compare the rates for horizontal events, where tau neutrino pileup is small. For example. 
Figures 22 (a-f ) show clear distinction between oscillation and no-oscillation scenarios, even 
in directions near horizontal, where there is no pile up. Furthermore, for E'"^ flux, where 



the pileup is very small [|12|, the ratio of ratios R discussed above ranges from 2.5 to 2.8. 
Thus, even without the tau neutrino pileup, the oscillation scenario can be distinguished 
from the no-oscillation scenario. 

The detection of z/^ oscillations with a point source might also be possible. With the 
resolution for the planned neutrino telescopes of 2°, the atmospheric background is reduced 
by 3.8 X 10^'^. For upward showers, this gives less than 1 event per year for -E^™ = 1 TeV, 
and even less for higher energy thresholds. Thus, if the point source has a flat spectrum, 
Fj^+p = 10~^^i?~^, then one would be able to detect tau neutrinos by measuring upward 
showers with E^^ = 1 TeV. In the more realistic case, when the point source has a steeper 
spectrum such as Sgr A* 0], a normalization of 10~''/cm^/s/sr/GeV would be 



sufficient for the detection of tau neutrinos with threshold of 1 TeV. Time correlations with 
variable point sources would further enhance the signal relative to the background. 

We have demonstrated that extragalactic sources of neutrinos can be used as a very-long 
baseline experiment, providing a source of tau neutrinos and opening up a new frontier in 
studying neutrinos oscillations. 
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APPENDIX A: TAU DECAY DISTRIBUTION 

The decay distribution of the tau neutrinos from tau decay has the following form, in 
terms of z = E^/Et-: 

fjji , 

-^ = T.M9o + P9\) ■ (Al) 

The polarization of the decaying r~ is P, which for neutrino V-A production of r~ is 
P = —1. The branching fraction into decay channel i is indicated by Bi. The distribution 
is normalized such that 

In Table 1, we show the functions go and gi for each decay mode, written in terms of z 
and Ti = mf/m^. Details of the calculational procedure can be found in Ref. p5| or in Ref. 



56| . For the multiprong tau decays, we approximate the distribution by a theta function, 



as indicated in the table. 



15 



REFERENCES 



[1] Super-Kamiokande Collaboration (Y. Fukuda et al.) Phys. Rev. Lett. 82, 2430 (1999); 
ibid, |hep-ex/9908(M)| ; Phys. Lett. B433, 9 (1998); Phys. Lett. B436, 33 (1998); ibid, 
Phys. Rev. Lett. 81, 1562 (1998). 

Kamiokande Collaboration, S. Hatakeyama et al, Phys. Rev. Lett. 81, 2016 (1998). 
Kamiokande Collaboration, K. S. Hirata et al, Phys. Lett.B205, 416 (1988); Phys. Lett. 
B280, 146 (1992); Y. Fukuda et al, Phys. Lett. B335, 237 (1994). 
D. Casper et al, Phys. Rev. Lett. 66, 2561 (1991); R. Becker-Szendy et al, Phys. Rev. 
D46, 3720 (1992). 

W. W. M. Allison et al, Phys. Lett. B391, 491 (1997). 
N. Fornengo, M. C. Gonzalez- Garcia and J. W. F. Valle, |hep-ph/0002T47 . 
MINOS Collaboration, E. Abies et al., "Main Injector Neutrino Oscillation Search", 
FERMILAB-PROPOSAL-P-875 (1995). 

K2K Collaboration, K. Nishikawa et al. "E362 KWK-PS proposal", March 1995; Nucl. 
Phys. (Proc. Supp.) B59, 289 (1997). 

ICARUS Collaboration, A. Rubbia et al, CERN-SPSLC-96-58 (1996). 
NOE Collaboration, M. Ambrosio et al, INFN-AE-98-09 (1998). 
OPERA Collaboration, K. Kodama et al, CERN/SPSC-98-25 (1998). 
S. Iyer, M. H. Reno, I. Sarcevic, Phys. Rev. D61, 053003 (2000). 

For a review of neutrino oscillations see, P. Fisher, B. Kayser and K. S. McFarland, 
Ann. Rev. Nucl Part. Scz. 49, 481 (1999). 
F. Halzen and D. Sahzberg, Phys. Rev. Lett. 81, 4305 (1998). 
T. J. Weiler, W. A. Simmons, S. Pakvasa and J. G. Learned, |hep-ph/9411432 . 



H. Athar, M. Jezabek and O. Yasuda, TMUP-HEL-0006 (2000); O. Yasuda, TMUP- 
HEL-0009 (2000). 

J. P. Rachen and P. Meszaros, Phys. Rev. D58, 123005 (1998). 

T. Stanev et al, |astro-ph/000"348^ . 

See, for example, URL |http : //antares . in2p3 . f r/. 



See, for example, URL |http : / /www . uao . gr^nestor/ 



See, for example, URL |http : / /amanda . berkeley . edu/ . 
See, for example, URL |http : / /amanda . berkeley . edu/km3/| . 

S. Bottai and F. Becattini, in Proceedings of the Cosmic Ray Conference (HE 4.2.15), 
Salt Lake City, Utah, August 17-23, (1999); S. Bottai and F. Becattini, |astro-| 
ph/0003179| . 



C. E. Fichtel, et al, Astrophys. J. Suppl 94, 551 (1994). 

M. Punch, et al. (Whipple Observatory Gamma Ray Collaboration) Nature (London) 
160, 477 (1992); A.D. Kerrick et al, Astrophys. J. 438, L59 (1995); J. Quinn, et al 
(Whipple Observatory Gamma Ray Collaboration) lAU Circular 6169 (June 16, 1995); 
J. A. Gaidos, et al. Nature (London) 383, 319 (1996); J. Quinn, et al, Ap. J. Lett. 
456, L63 (1996); J. H. Buckley, et al, Ap. J. 472, L9 (1996). 
[26] M. Catanese, et al, Ap. J. 501, 616 (1998). 

[27] For a discussion of the debate over the origin of high energy photons, see, e.g., J. H. 
Buckley, Science 279, 676 (1998); K. Mannheim, Science 279, 684 (1998); J. Rachen, 
|astro-ph/0003282 . 



[28] F. W. Stecker and M. H. Salamon, Space Sci. Rev. 75, 341 (1995). 



16 



[29] K. Mannheim, Astropart. Phys. 3, 295 (1995). 

[30] R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Astropart. Phys. 5, 81 (1996). 
[31] R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Phys. Rev. D58, 093009 (1998). 
[32] E. Waxman and J. N. Bahcall, |hep-ph/9909"286| ; C. D. Dermer, |astro-ph/0005440| . 
[33] R Meszaros and M. J. Rees, M. N. R. A. S. 258, 41P (1992); M. J. Rees and R Meszaros 
Ap. J. Lett. 430, L93 (1994). For a review, see T. Piran, Phys. Rept. 314, 575 (1999). 
[34] E. Waxman and J. N. Bahcall, Phys. Rev. Lett. 78, 2292 (1997). 

[35] E. Waxman and J. N. Bahcall, Phys. Rev. D59, 023002 (1999); J. N. Bahcall and E. 

Waxman, |hep-ph / 9902383| ; K. Mannheim, R. J. Protheroe and J. P. Rachen, |astro- 
I ph/9812398| . 

[36] G. Sigl, S. Lee, D. N. Schramm, and P. Coppi, Phys. Lett. B392, 129 (1997). 

[37] U. F. Wichoski, J. H. Macgibbon, and R. H. Brandenberger, 'High Energy Neutrinos, 



Photons and Cosmic Rays from Non-Scaling Cosmic Strings,' BROWN-HET-1115, [hep- 
ph/98054l9|. 



[38] R. J. Protheroe and T. Stanev, Phys. Rev. Lett. 77, 3708 (1996). 
[39] W. Rhode et al, Astropart. Phys. 4, 217 (1996). 

[40] R. Baltrusaitis et al, Astrophys. J. Lett. 281, L9 (1984); R. Bahrusaitis Phys. Rev. 
D31, 2192 (1985). 

[41] V. Agrawal, T. K. Gaisser, P. Lipari and T. Stanev, Phys. Rev. D53, 1314 (1996). 
[42] C. T. Hill and D. N. Schramm, Phys. Rev. D31, 564 (1985); S. Yoshida and M. Teshima, 

Prog. Theor. Phys. 89, 833 (1993); R. J. Protheroe and P. A. Johnson, Astropart. Phys. 

4, 253 (1996), Errat. ibid 5, 215 (1996). For a recent calculation of ultrahigh energy 

cosmic ray propagation, see, e.g., T. Stanev et al, Ref. [18]. 
[43] D. Seckel, T. Stanev and T. Gaisser, Astrophys. J 383, 652 (1991); G. Ingelman and M. 

Thunman, Phys. Rev. D54, 4385 (1996); C. Hettlage, K. Mannheim and J. G. Learned, 

AstropaH. Phys. 13, 45 (2000). 
[44] V. A. Naumov and L. Perrone, Astropart. Phys. 10, 239 (1999). 

[45] J. G. Learned and S. Pakvasa, Astropart. Phys. 3, 267 (1995); H. Athar, G. Parente and 
E. Zas, |hep-ph/0006123 ; J. Alvarez-Muniz and F. Halzen, "Detection of Tau Neutrinos 



in IceCube," unpublished. See also A. Husain, Nucl. Phys. Proc. Suppl. 87, 442 (2000) 

and |hep-ph/0004083| . 
[46] See, e.g., A. Breakstone et al., ISR collaboration, Phys. Rev. D30, 508 (1984). 
[47] The IceCube proposal can be found at |http : //pheno .physics .wise . edu/icecub4 
[48] A. Dziewonski, Earth Structure, Global, in Encyclopedia of Solid Earth Geophysics, ed. 

David E. James (Van Nostrand Reinhold, New York, 1989) pg. 331. 
[49] P. Lipari and T. Stanev, Phys. Rev. D44, 3543 (1991). 
[50] H. L. Lai et al. (CTEQ Collaboration), Eur. Phys. J. C12, 375 (2000). 
[51] L. Pasquah and M. H. Reno, Phys. Rev. D59, 093003 (1998). 

[52] The energy spectra of neutrino induced showers for E < 1 TeV, in the context of 
neutrino oscillations of atmospheric neutrinos, appears in T. Stanev, Phys. Rev. Lett. 
83, 5427 (1999). 

[53] J. G. Learned, private communications. 

[54] S. Markoff, F. Meha and I. Sarcevic, Astrophys. J. 522, 870 (1999). 
[55] T. K. Gaisser, Cosmic Rays and Particle Physics, (Cambridge University Press, Cam- 
bridge, 1990). 



17 



[56] P. Lipari, Astropart. Phys. 1, 195 (1993). 



18 



TABLES 



TABLE I. Functions go and gi in the tau neutrino energy {Ei,) distribution from r decays, in 
terms oi z = E^/Et- Note, X indicates hadrons, with X 7^ tt, p, a\. 



Process 



90 



91 



T — ^ i^T-ez/g 

T — ^ I^T-TT 
T VrP 

T — ^ Vt-O-I 



0.18 

0.12 
0.26 

0.13 
0.13 



I - 3.2 + |.3 



1 - 3.2 + §.3 
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FIGURES 




FIG. 1. Muon neutrino plus antineutrino fluxes for AGN models (blue lines, upper curve 
at low energy corresponds to AGN_M95, while the lower curve is for AGN_SS model), 
GRB (green line), topological defects models (black lines, upper curve corresponds to the 
TD_WMB, while the lower curve is for TD_SLSC), flux (lower red line at low energy) 
and E'"^ (upper red line at low energy) and angle-dependent atmospheric (ATM) flux (green 
shaded area). 
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E.(GeV) E,(GeV) 

FIG. 2. Muon neutrino plus antineutrino flux (black line), the effect of its attenuation for 
= 0° (red line) and tau neutrino plus antineutrino upward flux for the same initial flux and 
the same nadir angle (blue line) for a) E'^ flux b) flux c) AGN_SS and d) AGN_M95. 
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E.(GeV) E,(GeV) 

FIG. 3. Muon neutrino plus antineutrino flux (black line), the effect of its attenuation 
for 6 = 0° (red line) and tau neutrino plus antineutrino upward flux for the same initial flux 
and the same nadir angle (blue line) for a) TD_WMB, b) TD_SLSC, c) GRB_WB and d) 
the atmospheric flux ATMOS. 
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FIG. 4. The energy dependence of the muon neutrino flux (red hne) and the tau neutrino 
flux (blue hne) for nadir angle 6 = 0°, 9 = 30° and = 60° normalized to the initial flux for the 
AGN_M95. 
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FIG. 5. Diagrams for neutrino interactions contributing to the muon production. 
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FIG. 6. Muon event rate as a function of nadir angle for energies 1 TeV, 10 TeV and 100 
TeV. Muon rates including the contribution from tau decay (blue line) compared with the 
background from muon neutrinos (red line) for a) flux and b) flux. 
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FIG. 7. Muon event rate as a function of nadir angle for energies 1 TeV, 10 TeV and 100 
TeV. Muon rates including the contribution from tau decay (blue line) compared with the 
background from muon neutrinos (red line) for a) AGN_SS and b) AGN_M95. 
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FIG. 8. Muon event rate as a function of nadir angle for energies 1 TeV, 10 TeV and 
100 TeV. Muon rates including contribution from tau decay (blue line) compared with the 
background from muon neutrinos (red line) for a) TD_WMB and b) TD_SLSC. 
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FIG. 9. Muon event rate as a function of nadir angle for energies 1 TeV, 10 TeV and 100 
TeV. Muon rates including the contribution from tau decay (blue line) compared with the 
background from muon neutrinos (red line) for GRB_WB. 
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FIG. 10. Muon event rate as a function of nadir angle for energies 1 TeV, 10 TeV and 
100 TeV for a) Atmospheric muon neutrino and b) Atmospheric tau neutrino. 
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FIG. 11. Diagrams for neutrino interactions contributing to the shower events. 
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FIG. 12. Hadronic/EM event rates as a function of nadir angle for =1 TeV, 10 TeV 
and 100 TeV. Hadronic/EM event rates from (blue line) compared hadronic/EM event 
rates from z/^ plus z/g (red line) for a) and b) E""^. 
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FIG. 13. Hadronic/EM event rates as a function of nadir angle for =1 TeV, 10 TeV 
and 100 TeV. Hadronic/EM event rate from z/^ (blue line) compared hadronic/EM event 
rate from z/^ plus Ue (red line) for a) AGN_SS and b) AGN_M95. 
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FIG. 14. Hadronic/EM event rates as a function of nadir angle for = 1 TeV, 10 TeV 
and 100 TeV. Hadronic/EM event rates from Uj. (blue line) compared hadronic/EM event 
rates from plus (red line) for a) TD_WMB and b) TD_SLSC. 
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FIG. 15. Hadronic/EM event rates as a function of nadir angle for = 1 TeV, 10 TeV 
and 100 TeV. Hadronic/EM event rates from z/,- (blue line) compared hadronic/EM event 
rates from z/^ plus z/g (red line) for GRB_WB. 
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FIG. 16. Hadronic/EM event rates as a function of nadir angle for = 1 TeV, 10 TeV 
and 100 TeV. Hadronic/EM event rates from 1/^ (blue line) compared hadronic/EM event 
rates from i/^ plus z/g (red line) for a) atmospheric muon neutrino and b) atmospheric tau 
neutrino. 
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FIG. 17. Ratio of Hadronic/EM event rate of z/^ plus plus v^, assuming oscillation 
scenario and v^^ plus z/g in the standard model as a function of nadir angle for energies 1 
TeV, 10 TeV and 100 TeV for a) E"^ and b) E'"^ flux. 
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FIG. 18. Ratio of Hadronic/EM event rates of plus plus z/g assuming oscillation 
scenario and z/^ plus z/g in the standard model as a function of nadir angle for energies 1 
TeV, 10 TeV and 100 TeV for a) AGN_SS and b) AGN_M95 flux. 
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FIG. 19. Ratio of Hadronic/EM event rate of v^. plus z/^ plus v^. assuming oscillation 
scenario and z/^ plus z/g in the standard model as a function of nadir angle for =1 TeV, 
10 TeV and 100 TeV for a) TD_WMB and b) TD_SLSC flux. 
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FIG. 20. Ratio of Hadronic/EM event rate of Vt plus v^j, plus z/g assuming oscillation 
scenario and plus z/g in the standard model as a function of nadir angle for energies 1 
TeV, 10 TeV and 100 TeV for GRB_WB. 
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FIG. 21. Ratio of Hadronic/EM event rate of Vr plus plus z/g assuming oscillation 
scenario and plus in the standard model relative to the ratio of the muon event rate 
of Vr plus assuming oscillation scenario and z/^ in the standard model (see Eq. (15)) as 
a function of nadir angle for threshold energies 1 TeV, 10 TeV and 100 TeV. 
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FIG. 22. Ratio of Hadronic/EM event rate to muon event rate for the oscillation and 
no-oscillation scenarios as a function of nadir angle for threshold energies (a-b) 1 TeV, (c-d) 
10 TeV and (e-f) 100 TeV for the indicated fluxes. 
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